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ABSTRACT 

GRS 1915+105 is a widely studied black hole binary, well known because of its ex- 
tremely fast and complex variability. Flaring periods of high variability alternate with 
"stable" phases (the plateaux) when the flux is low, the spectra are hard and the tim- 
ing properties of the source are similar to those of a number of black hole candidates 
in hard spectral state. In the plateaux the power density spectra are dominated by 
a low frequency quasi periodic oscillation (LFQPO) superposed onto a band limited 
noise continuum and accompanied by at least one harmonic. In this paper we focus on 
three plateaux^ presenting the analysis of the power density spectra and in particular 
of the LFQPO and its harmonics. While plotting the LFQPO and all the harmonics 
together on a frequency-width plane, we found the presence of a positive trend of 
broadening when the frequency increases. This trend can shed light in the nature of 
the harmonic content of the LFQPO and challenges the usual interpretation of these 
timing features. 



1 INTRODUCTION 



X-ray binaries are binary systems where a stellar compact 
object, a neutron star (NS) or a black hole (BH) accretes 
matter from a companion star, producing bright X-ray emis- 
sion. Since its launch in 1995, the Rossi X-Ray Timing Ex- 
plorer (RXTE) has allowed a direct comparison between the 
spectral and timing properties of large sample of XRBs, 
which allowed the definition of the typical phenomenolog- 



ical properties of sources hosting a NS or a BH (see Psaltis 
12006 for a review). On the base of those properties, quite a 



number of sources have been identified as possibly hosting a 
BH (black-hole candidates, BHC): for a number of systems 
this prediction has been confirmed via dynamical studies. 
Most BHCs have a low-mass companion star and are tran- 
sient sources (black-hole transients, BHT) showing an evo- 
lution of their spectral/timing properties during their out- 
bursts. This evolution has been described by many authors 
as a transition through a combination of 3-4 states (see 
Remillard fc McClintock|2006| and |Belloni|2010] for a review) 



mainly defined by the shape and features in the power den- 
sity spectra (PDS) and by the hardness of the energy spec- 
tra. The typical PDS are composed by a band-limited noise 
continuum that varies during the outburst and by a number 
of quasi-periodic oscillations (QPOs) observed in association 
with certain states both at low (mHz to few tens of Hz) and 
high (tents to hundreds Hz) frequencies ( [van der KlispOOG] ). 
Various models have been proposed in order to explain the 
quasi periodic variability and the noise continuum (see 



|der Klis]|2005| for a review, [Ingram fc Done|2011| for a more 
recent model) but none of them, at the present time, can 



give a comprehensive explanation of all the observed phe- 
nomena. 

In this paper we focus on the low frequency QPO (LFQPO) 
of GRS 1915-1-105. LFQPOs are common features in the 
PDS of BHTs (unlike high-frequency QPOs, with frequen- 
cies larger than ~ 30 Hz, that are rather rare) during their 
hard spectral states. They appear as strong peaks at a 
centroid frequency between 0.01 and 15 Hz, with a frac- 
tional rms amplitude reaching values >15% in sources like 
XTE J 1550- 564 and GRS 1915-^105 ( |Sobczaket al.|2000[ 
Reig et al. 20001. LFQPOs are often accompanied by one 



or more further QPO peaks whose frequency is in harmonic 
relation with that of the LFQPO itself. Those are tradition- 
ally interpreted as resulting from the Fourier decomposi- 
tion of the quasi-periodic signal responsible for the LFQPO. 
Nonetheless, questions about the identification of the fun- 
damental peak and the genuineness of the harmonic rela- 
tionship of the peaks have been raised in the recent works 
of |Raoet al.jpMot and [Rodriguez fc Varnierel ( [20TT| . With 
the aim of investigating the relation between the LFQPO 
and its harmonics, we analysed a sample of RXTE obser- 
vations from GRS 1915+105, which is particularly favorable 
for this study due to the strength of the main QPO peak and 
the number of the harmonics. The source is a largely studied 
XR B (see [Fender fc Be"IIom|[2004| for a review) with a K-M 
HI ( Greiner et al.|2001b I star donating matter to a dynami- 
cally confirmed BH ( Greiner et al.|2001a \ . It was discovered 
as a transient source on August 15, 1992 by the WATCH m- 



strument on board on the GRAN AT saielWte [ Castro- Tirado 
et al.|1992| and since then it always remained bright in the 
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Table 1. List of the RXTE/PCA archival observations analysed 
in this paper. The observations are divided in three blocks, cor- 
responding to the three subsequent plateaux in Fig. IT] 



Obs 



MJD 
(50000.0) 



Counts 
(Crab units) 



10408- 


-01- 


-22- 


-01 


275.2230 


0.818 


10408- 


-01- 


-22- 


-02 


275.3564 


0.796 


10408- 


-01- 


-23- 


-00 


278.4915 


0.873 


10408- 


-01- 


-24- 


-00 


280.1696 


0.790 


10408- 


-01- 


-25- 


-00 


283.4941 


0.732 


10408- 


-01- 


-27- 


-00 


290.5768 


0.715 


10408- 


-01- 


-28- 


-00 


298.5349 


0.699 


10408- 


-01- 


-29- 


-00 


305.3730 


0.709 


10408- 


-01- 


-30- 


-00 


313.3129 


1.141 


10408- 


-01- 


-31- 


-00 


320.1960 


0.922 


20402- 


-01- 


-04- 


-00 


415.1299 


0.864 


20402- 


-01- 


-05- 


-00 


421.9758 


0.551 


20402- 


-01- 


-07- 


-00 


436.6608 


0.511 


20402- 


-01- 


-08- 


-00 


441.9193 


0.519 


20402- 


-01- 


-08- 


-01 


442.1189 


0.485 


20402- 


-01- 


-09- 


-00 


448.2834 


0.427 


20402- 


-01- 


-10- 


-00 


455.9928 


0.387 


20402- 


-01- 


-11- 


-00 


462.0617 


0.356 


20402- 


-01- 


-12- 


-00 


471.0688 


0.346 


20402- 


-01- 


-13- 


-00 


477.8727 


0.368 


20402- 


-01- 


-14- 


-00 


480.8806 


0.357 


20402- 


-01- 


-15- 


-00 


488.7786 


0.317 


20402- 


-01- 


-16- 


-00 


501.8847 


0.312 


20402- 


-01- 


-18- 


-00 


512.8887 


0.329 


20402- 


-01- 


-19- 


-00 


517.0462 


0.294 


20402- 


-01- 


-20- 


-00 


524.9231 


0.315 


20402- 


-01- 


-21- 


-00 


533.8351 


0.309 


20402- 


-01- 


-49- 


-00 


729.3267 


0.759 


20402- 


-01- 


-49- 


-01 


730.3943 


0.726 


20402- 


-01- 


-50- 


-01 


737.4041 


0.600 


20402- 


-01- 


-51- 


-00 


743.2941 


0.597 


20402- 


-01- 


-52- 


-00 


746.5504 


0.602 



X-ray sky (~ 0.5 - 2 Crab), undergoing an 18-year long out- 
burst that at the time of writing is stiU ongoing. Starting 
from 1996, GRS 1915-f 105 has been monitored with RXTE. 
The hght curve in Fig. [l] is an example of the long-scale 
behavior of the source: quiet periods of low flux and low 



variability lasting hours to tens of days (the plateaux, Fender 
|et al.|1999l ) , alternate with flaring activity phases of extreme 
variability, when all the X-ray properties, flux, spectral pa- 
rameters and PDS features can vary on timescales down 
to the ms ( [Greiner et al.||1996| |Bellom et aL]|1997[ |Belloni| 
et al. [[2000 1. While peculiar and complex during the flaring 
phases, the behaviour of the source during the plateaux is 
similar to that of many other BHTs. In particular, the PDS 
show the typical band limited noise complex (extending up 
to ~ 100 Hz) and a strong LFQPO with several harmonics 
| |Trudolyubo"v|[200H [Reig et al.[[2000| We analysed RXTE 
data from the three plateaux in Fig. [ll the data sample and 
the reduction method are described in Sect. [2J while the 
modeling of the PDS is described in Sect. [3] Sect.[4]presents 
our results, which are discussed in the final Sect.js] 




100 200 300 +00 500 600 700 
Time (MJD-50000.0) 

Figure 1. The RXTE/PCA observations analysed in this paper 
(black points) superimposed onto the RXTE/ ASNL daily— average 
light curve between October 1995 and December 1997 (from [Bel^ 
jloni et al.j ( |2000| ). The flux level is different because of the dif- 
ferent energy band of PCA and the ASM. PCA data belong to 
three plateaux, the middle one being significantly longer that the 
others. 




Frequency (Hz) 

Figure 2. The most general multi-Lorentzian model we applied 
to the PDS. The Lorentzian components are labeled according to 
[Belloni et al.| l [2002| . The hjy component accounts for residuals 
in the fit in a few observations, but is never significant. 



2 OBSERVATIONS AND DATA REDUCTION 

We analysed a set of thirty-two RXTE/PCA observations 
from the RXTE NASA's archive (Table [I]) performed be- 
tween July 1996 and October 1997 during the three subse- 
quent plateaux shown in Fig, [l] ( for a m ore general studies 
of the same plateaux, Trudolyubov|2001 1 . 
Power Density Spectra (PDS) have been extracted in a 
~ 2 — 14 keV energy band from PCA light curves (single 
bit data with temporal resolution 125/^s) with a Nyquist 
frequency of 2048 Hz, normalised according to [Leahy et al.[ 
( 1983 1 and converted to square fractional root mean square 
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Frequency (Hz) 




Frequency {Hz) 



Figure 3. PDS from observation 20402-01-19-00 (A) and 10408-01-27-00(B) with our best fitting model. Continuum components 
are plotted as dashed lines, QPOs with filled lines. The high frequency component L; (see Fig.[2| that is only observed during the long 
plateau ||Trudolyubov|2001[l can be seen in panel (A). 



(rms) deviation ( Belloni fc Hasinger||1990 l. In order to im- 
prove the statistics, we obtained a PDS for each observation 
as the average of the Fourier spectra extracted from con- 
secutive 16 — s segments of the light curve. The averaged 
PDS have been re-binned logarithmically in frequency and 
fitted with the standard Xspec (version 11.3.2) fitting pack- 
age using a diagonal response, therefore performing a simple 
direct y^ minimization. We adopted a model composed by 
Lorentzian functions only (see Fig. pi following |Belloni et al 
(20021. As argued by the authors, this approach has the ad- 



vantage that the spectrum is described by components that 
are directly comparable one to the other, with no assump- 
tion about the origin of each one. A component has been 
considered significant if N/gm > 3, where A'' is the normali- 
sation and (Tjv its uncertainty, obtained with the chi-squared 
minimization. The errors on individual powers of the fitted 



PDS were computed following van der Klis (19891. We ac- 



counted for the counting statistic in PDS using an automatic 
subtraction of Poissonian background, based on the estimate 
of PCA dead time by Zhang et al. ( 1995 1. A flat component 



extending over all the frequencies appears in a few PDS as a 
result of a non-perfect Poissonian subtraction. We corrected 
for this effect while fitting the PDS. 



3 MULTI-LORENTZIAN MODELING OF THE 
PDS 

Figure |3] shows two examples from the sample of PDS we 
analysed and the associated best-fitting model. Depending 
on the observation, five to eight Lorentzians are required 
for a good fit ( x^ < 1-4). We interpret each Lorentzian as 
a different spectral component (with the exception of dou- 
ble peaked QPOs, see below) assigning labels as in Fig. [2] 
We define the frequency u and width A of a component 
as the central frequency and the FWHM of the associated 
Lorentzian function. 
One to three Lorentzian components account for the band 



limited continuum (see Fig. [2] and [3| : a low frequency one, 
Li,, appears in 29 over 32 observations at j/f, < 0.2 Hz. Lh is 
often required at a centroid frequency Vh , approximately lo- 
cated under the LFQPO Llf- A high frequency component 
Li only appears in observations from the longest plateau in 
Fig. [l] at z/j > 42 Hz. For a study of this component and 
its connection with "short" and "long" plateaux, see |Tru-| 
|dolyubov|[200T| 

One Lorentzian is usually enough to fit the LFQPO peak, 
unless it happens to be double-peaked. The presence of a 
double peak is possible because the QPO is known to drift 
in frequency on time scales shorted than the average obser- 
vation (see e.g. Markwardt et al. 1999 ). When this is the 



case, we fit two Lorentzians to the QPO and we take the 
sum-function as a single spectral component Llf to rep- 
resent the LFQPO. We consider the frequency of the most 
significant (larger %rms) peak as the frequency oi Llf, vlf 
and define its width A_lf as the FWHM of the sum-function 
of the overlapping Lorentzians. The error on Alf is analyt- 
ically computed from the equation of the sum-function. The 
quality factor Q — v/A (typical indicator for the coherence 
of a signal, a QPO is traditionally defined by Q >2) of Llf 
is always above 3, reaching a maximum value of 10.4 (see 
Tabled). 

Although this is not the main topic of this paper, it is worth 
to mention that our data are consistent with the relations 
between Vh, vi and ulf found in Belloni et al. (20021 for a 



sample of sources not including GRS 1915+105. 
After fitting for the main QPO peak and the broad band- 
limited noise components Lb, Lh and Li, further peaks or 
bumps in the PDS are best fitted by up to three compo- 
nents in harmonic relation with Llf, labeled L//, Lm and 
Li/2- Their central frequencies are vu ~ 2vlf, viii ~ 'ivLF 
and Vi/2 ~ O.^vlf respectively. In three observations a fur- 
ther component {Liv in Figure [2| needs to be added to 
our model in order to take into account residuals in the fit, 
whose centroid frequency is consistent with ~ ^vlf- How- 
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a measure for the break frequency of a broad Lorentzian and 
around this frequency the component contri butes most of 
its power per logarithmic frequency intervarj ([Belloni et al 



Figure 4. Top panel: Central frequency oi Li^p (black stars) , Ljj 
(empty stars), Lnj (black circles), L^ /2 (black diamonds) and Lh 
(empty diamonds) versus the frequency of Lj^p. The black circles 
in an empty frame are 2 a detections of Ljij. Lines of constant 
ratio R= 0.5, 1, 1.5, 2 and 3 are plotted: the solid line is R=l. 
Bottom panel: same as top panel, but the characteristic frequency 
i^max = yji^^ + (A/2)-^ is plotted instead of u. 



ever, this component never results significant (slightly less 
than 2(T detection) and therefore we did not consider it in 
our subsequent analysis. The component Lu always appear 
as a rather coherent peak, with 2 < Q < 10. Lm and I/1/2 
instead are often broad, with Q ranging from less than 1 up 
to ~8 (see Table Isl III [^ . 



3.1 Harmonic relations 

Figure [4] (top panel) shows the frequencies vu, um, vlf, 
fj/i and Uh plotted against the corresponding vlf (the data 
are reported in the tables from [2] to [m . The harmonic rela- 
tion between the Llf and Lu is evident from the figure. 
Quite some scattering is observed instead in the case of the 
broader components Lm and I/1/2, with a tendency of the 
points to lie slightly below the harmonic relation line. Fig- 
ure 14] (bottom panel) suggests that the harmonic ratio is 
better represented for Lm and I/1/2 if we consider, instead 
of the frequency of the compone nts, their 'character istic fre- 
quency' fmax = ^Jv"^ + (A/2)2 ( [Belloni et al.|2002| . This is 



max of the com- 
by the authors) is in a 3/2 ratio 



|2002|. When considering v (upper panel of Fig. |4|| , a fit to 
the Lm and I/1/2 data-points with a line of constant ratio 
V ^ AvLF gives A = 2.91 ± 0.02 (x^= 131.2, 24 d.o.f) and 
A = 0.455 ± 0.004 (x^= 412.3, 27 d.o.f) for the two compo- 
nents respectively. The fitted slope is inconsistent with the 
expected harmonic ratios between the frequency of the each 
component and that of the Llf ■ If we consider Vmax instead 
(Fig. [i] bottom panel), the fit gives A = 3.00 ± 0.02 for Lm 
{X^= 116.4, 24 d.o.f) and A = 0.49 ± 0.004 for L1/2 (x^= 
191.4, 27 d.o.f). In this case, the fitted slope is consistent 
with the expected harmonic ratios on the la level for Lm 
and on the 3a level for I/1/2. Although the fits are poor due 
to the scatter between the data-points, they support the vi- 
sual impression given by Fig. |4] that the harmonic ratio of 
the Lm and I/1/2 is better represented when Vmax is used 
instead of v. 

The frequency of Lh is included in Figure [4] following a re- 
cent s tudy of the PDS of XTE J1550-564 performed by [Rao] 
et al.| (2010) . The authors found that the u, 
ponent Lh (named L- 
with h'max(LF), being at three times the characteristic fre- 
quency of the sub-harmonic I/1/2 ■ We find that, in the case of 
GRS 1915-1-105, Vmaxih) lies slightly below the line indicat- 
ing a constant ratio R=3/2 with VmaxiLF) (Figure |4| bottom 
panel). As Lh is a rather broad component, this discrepancy 
is even stronger if we consider the centroid frequencies Vh 
and ulf (Figure |4] top panel). 



4 HARMONICS IN A WIDTH-FREQUENCY 
PLANE 

Fig. [5] shows the width of the Llf and its higher order har- 
monics Lu and Lm plotted against the centroid frequency 
of the Llf- The three components describe three clear tracks 
on the diagram, and the width of each component grows 
with frequency. Moreover, it seems that the highest is the 
order of the harmonic, the broader the peak is. Nonetheless, 
a different interpretation arises when considering the plot 
of Fig. [6] which provides a deeper insight on the behaviour 
of the harmonics. Here Llf, Lu and Lm are considered 
again on a width-frequency plane, but the width A of each 
component is plotted against its own frequency v, without 
implying a relation with the Llf (the data are reported in 
the tables from|2]to|6|. Lines of constant quality factor Q 
are also plotted. Different symbols distinguish harmonics of 
different orders, with the same key used in Fig. |5] With the 
exception of a few outliers, all points follow the same trend 
on this plot, where coherence decreases as centroid frequency 
increases. The harmonics can not be distinguished based on 
their width, as the tracks from different harmonics overlap 
one with another. Rather than depend on the order of the 
harmonic, a given A seems to be associated with a given 
frequency. A similar trend is observed if Vmax is plotted in- 
stead of v. 



Note that u„ 



becomes equal to v for narrow components. 
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Figure 5. For each L^p (black stars), Ljj (empty stars) and 
Ljij (black dots) detection, the width is plotted against the cen- 
troid frequency of the Li^p in the same PDS. The black dots in a 
white frame are 2 a detections consistent with Ljjj. The dashed 
line shows the slope of a constant quality factor (Q=10) trend. 



The relation between the quality factors of the different har- 
monics of a non-sinusoidal signal is determined by the nature 
of the quasi-periodicity of the signal. The width-frequency 
plots in Fig.lTlhave been produced by assuming the u of each 
harmonic peak in our data and calculating the expected A in 
case of a purely frequency - or purely amplitude - modulated 
signal. For a purely amplitude modulation, the A of all the 
harmonic peaks would be the same as the amplitude of the 
Llf- On a frequency- width diagram that would result in a 
series of parallel tracks, one per each harmonic as shown in 
Fig. [7] panel (a). The diagram is clearly different than that 
resulting from the real data in Fig. |6] In the case of a fre- 
quency modulation, all the peaks have the same quality fac- 
tor, that of the Llf- This case is shown in Fig. JTlpanel (b), 
which again is not consistent with the real data: although a 
single track is formed, it remains parallel to the Q — const 
lines, not reproducing the curved shape that appears in the 
real GRS 1915-1-105 data set. A width-frequency diagram 
consistent with a pure frequency modulation was found in- 
stead for XTE J1550-564 ( |Rao et al.|2"010l ). 
The frequency-width diagram for components I/1/2 and Lh 
is shown in Fig. [8] A comparison of Fig. |8] and |6] shows that 
Li/2 and Lh do not share the behaviour of the other har- 
monics. 



5 DISCUSSION AND CONCLUSIONS 

We have analysed the PDS of GRS 1915-f 105 during three 
plateaux states. The overall shape of the PDS is that typical 
of BHTs in their hard-intermediate state, with a strong low 
frequency QPO peak Llf superimposed to a band limited 
noise continuum (Fig. [3][2|. The spectra can be fitted with 
a combination of Lorentzian components, several of which 
are in harmonic relation with Llf- We detected up to two 
higher order harmonics (L// and Lm) forming an harmonic 
series of 1:2:3 with Llf, plus the components I/1/2 and Lh 
with frequency close to ~ 0.5 and ~ 1.5 times the frequency 
of the Llf respectively. Few observations show residuals at 



frequency consistent to ~ 4:Llf, but a further Lorentzian 
component added at this frequency never results significant 
above 2a. 

In the case of XTE J 1550-564, |Rao et al.| ( |2010| ) proposed 
Li/2 as the fundamental frequency of an harmonic 1:2:3:4 
series including I/1/2, Llf, Lh and Lu- In a similar scenario 
for GRS 1915-f 105, Lm should also be included, leading to 
a 1:2:3:4:6 series. Nonetheless, there is no strong evidence of 
a 2/3 harmonic ratio between the frequency of the Llf and 
of Lh in the case of GRS 1915-f 105, even when the char- 
acteristic frequency of the components Umax is considered 
in place of v. Moreover, the plot in Figure p] (which is not 
significantly affected by the choice of f or fmax) indicates 
that Llf and its higher order harmonics are related beyond 
their frequency ratio as they broaden and lose coherence 
together as their frequency increases, while I/1/2 and Lh be- 
have differently (Fig. [8| as if they did not belong to the Llf 
harmonic series. We conclude that the I/1/2 is likely not the 
fundamental in the Llf harmonic series, i.e. is not included 
in the series itself. Z/1/2 and Lh could instead be produced by 
a different phenomenon than that responsible for the Llf- 
Although the plot in Fig. |6] evidences a relation between 
Llf and its higher order harmonics Ljj and Ljjj, it also 
rises doubts on the nature of this relation. The commonly 
accepted idea is that the Llf and the harmonics describe 
together the same signal from a quasi-periodic oscillator, ap- 
pearing as separated peaks in the PDS as a consequence of 
the Fourier representation. However, it is not trivial to iden- 
tify a signal modulation that is able to produce the trend 
in Fig. [6] For example, we have shown that such a behavior 
cannot be reproduced by a simple combination of frequency 
or amplitude modulation. Moreover, Fig. [4] shows that the 
frequency of Lm across different observations is better dis- 
tributed around an harmonic ratio of 3 with the frequency of 
Llf when considering lymax instead of i^Fl This is not what 
is usually expected in the context of harmonic decomposi- 
tion, where the centroid v is the frequency expected to be 
in harmonic ratio with the fundamental frequency. A broad 
component whose Umax is in harmonic ratio with the funda- 
mental frequency can not be regarded as a proper harmonic, 
but only as a signal whose characteristic frequency is con- 
sistent with an harmonic of the fundamental frequency. The 
complexity of these results invites to consider a scenario in- 
volving more than one real oscillator, where a physical phe- 
nomenon is triggering oscillations at multiple frequencies in 
the accretion disk or in the corona, each one resulting in 
an harmonic component of the PDS. The trigger could be 
such that the life time (i.e. the coherence) of each oscilla- 
tor is determined only by its own frequency, in agreement 
with the the trend in Fig. [6] This toy-model example gives 
a feeling of the new perspectives that would be opened in 
the interpretation of the PDS if the harmonics were proven 
to have a meaning beyond the Fourier representation of a 
single quasi-periodic signal. Either way, the interpretation 
of the width-frequency plane of GRS 1915-f 105 promises to 
offer new insights on the QPO generation phenomenon in 
this system and in other XRBs. 



^ Note that, because Llh is narrow, there is no significant dif- 
ference between Vmax and f for this component. The same holds 
for Li[. 
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Figure 6. Harmonics on a width-frequency plane. Symbols are the same as in Figurelm for each L^p (black stars), Ljj (empty stars) 
and Ljij (black dots) component detected in the PDS, the width is plotted against the centroid frequency. The black dots in a white 
frame are 2 a detections consistent with Ljij. Lines of constant quality factor Q = I'/A are plotted (dashed lines, Q increasing from 1 
to 10 going from top to bottom). 
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(a) (b) 

Figure 7. Panels (a) and (b) show how the plot in Fig. [6] would look like in case the harmonics origin from a signal modulated in 
amplitude and frequency respectively. The points have the same frequency as the real data, but the width of each harmonic is : (a) the 
same as the width of the Llf in the same observation (amplitude modulation) (b) calculated in order to obtain a constant Q for all the 
harmonics from the same PDS (frequency modulation). 
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Obs ID 



^LLF 



LogCi-CHz)) 



Figure 8. Width versus centroid frequency (same as Fig.pl for 
tlie sub-liamionic Li /2 (black diamonds) and tlie component Lh 
(empty diamonds). Lines of constant quality factor Q = i//A are 
plotted as in Fig. ml (dashed lines, Q increasing from 1 to 10 from 
top to bottom) 
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10408-01-23-00 
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10408-01-25-00 


1.121 


10408-01-27-00 
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1.045 


20402-01-51-00 


1.389 


20402-01-52-00 
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-1-0.005 
-0.006 
+0.007 
-0.007 
+0.005 
-0.010 
+0.002 
-0.003 
+0.002 
-0.002 
+0.002 
-0.002 
+0.003 
-0.004 
+0.010 
-0.010 
+0.006 
-0.007 
+0.021 
-0.025 
+0.004 
-0.005 
+0.006 
-0.006 
+0.011 
-0.010 
+0.010 
-0.011 
+0.007 
-0.007 
+0.006 
-0.006 
+0.008 
-0.009 
+0.006 
-0.007 
+0.007 
-0.008 
+0.013 
-0.016 
+0.005 
-0.004 
+0.008 
-0.009 
+0.011 
-0.011 
+0.005 
-0.005 
+0.007 
-0.007 
+0.012 
-0.013 
+0.006 
-0.007 
+0.011 
-0.014 
+0.002 
-0.004 
+0.002 
-0.003 
+0.003 
-0.005 



359+''°0^ 

0.345irol^ 
0.468irol'8 

'-'■^ '"-0.024 
,-, 179+0.005 
>J.1(^_0.007 
(,(,-,+0.005 
U.U^1_0 004 

'-'■J^''^_0.003 
n ncri +0.009 

1 fiQ9+0.034 
1-DO^_0.034 
439+0.023 
u.*o»_g 023 

700+0081 
U.(UU_o 074 

0.466±rof9 

0.819irof3 

668+0035 
U.DDO_o 052 

n oor,+0.038 

u.oo^_0 039 

544+0029 
U.0t1_o 033 

625+0026 
U.DZO_o 024 

o.eoiirofs 

0.483irof2 

Q40+0023 
U.y4U_o 035 

749+0.054 
u./ty_o 067 

O QQO + 0.016 

u.o»o_o 016 

548+0036 
U.OtO_o 037 

489+0063 
U.y:Oa_o 054 

r. oqo + 0.021 
U.OOO_o 019 

485+0039 
U.400_o 037 

fi56+0109 
U.DOD_o 067 

291+0019 
u.zai_o 018 

310+0039 
U.01U_o 031 

154+0008 
U.101_o 010 

n 1fi«+0008 
U.16«_ooo9 
I3fj+0.010 
U.lciD_o,oi2 



Table 2. Frequency ullf, width Ai^^p 
QPO Llp from the best fit to the PDS. 
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ObsID 



>^ii 



Obs ID 



i^iii 



10408-01-22-01 
10408-01-22-02 
10408-01-23-00 
10408-01-24-00 
10408-01-25-00 
10408-01-27-00 
10408-01-28-00 
10408-01-29-00 
10408-01-30-00 
10408-01-31-00 
20402-01-04-00 
20402-01-05-00 
20402-01-07-00 
20402-01-08-00 
20402-01-08-01 
20402-01-09-00 
20402-01-10-00 
20402-01-11-00 
20402-01-12-00 
20402-01-13-00 
20402-01-14-00 
20402-01-15-00 
20402-01-16-00 
20402-01-18-00 
20402-01-19-00 
20402-01-20-00 
20402-01-21-00 
20402-01-49-00 
20402-01-49-01 
20402-01-50-01 
20402-01-51-00 
20402-01-52-00 



5.386 
5.024 
6.794 
4.399 
2.194 
1.250 
1.879 
3.780 
8.121 
8.186 
8.586 
5.582 
6.195 
7.784 
6.817 
5.590 
5.818 
5.733 
5.570 
7.242 
6.966 
4.541 
5.949 
6.437 
4.424 
6.403 
7.100 
5.178 
5.387 
2.049 
2.723 
2.771 



+0.051 
-0.057 
+0.036 
-0.038 
+0.079 
-0.106 
+0.025 
-0.035 
+0.009 
-0.010 
+0.005 
-0.006 
+0.008 
-0.008 
+0.020 
-0.021 
+0.154 
-0.206 
+0.066 
-0.067 
+0.286 
-0.358 
+0.019 
-0.010 
+0.039 
-0.037 
+0.051 
-0.043 
+0.048 
-0.047 
+0.039 
-0.031 
+0.022 
-0.023 
+0.041 
-0.038 
+0.026 
-0.026 
+0.031 
-0.030 
+0.046 
-0.049 
+0.015 
-0.016 
+0.046 
-0.057 
+0.056 
-0.053 
+0.015 
-0.015 
+0.034 
-0.035 
+0.061 
-0.065 
+0.040 
-0.043 
+0.035 
-0.038 
+0.011 
-0.010 
+0.007 
-0.008 
+0.012 
-0.013 



1 601+"-l^* 
-, 0-O+0.151 

l.U(0_g j^28 



1.110 



1+0.164 
-0.122 



0.561«:»i 



0.733" 



0.067 



*-OOO_0.408 

lU.4y^_g ggg 

l.iU»_o.o99 

2.012«;1« 

+0.286 
-0.260 
+0.146 



1.998 
1.585 



1 344+0.222 

1.0'i1_g QgQ 

1.702t°o:\ll 

1 4fi2+"-"5 

1.4DZ_g j^j^g 

2 626+°°'5'' 

Q4g+0.095 

1 006+°°™ 
1 454+0-253 

i.76it°:f,l 

0.808t°:°?? 

a+0.171 

-0.283 

1 QQc;+0.307 

±.yaO_g 3Qg 

.+0.159 
"'-0.132 
7+0.062 
'-0.100 
n •308+0.049 

n 492+0.040 

n 4SS+0-049 
U.488_o.043 



1.688] 
1.995^ 
0.544^ 
0.577"^ 



10408-01- 
10408-01- 
10408-01- 
10408-01- 
10408-01- 
10408-01- 
10408-01- 
10408-01- 
20402-01- 
20402-01- 
20402-01- 
20402-01- 
20402-01- 
20402-01- 
20402-01- 
20402-01- 
20402-01- 
20402-01- 
20402-01- 
20402-01- 
20402-01- 
20402-01- 
20402-01- 
20402-01- 
20402-01- 



22-01 

22-02 

23-00 

24-00 

25-00 

27-00 

28-00 

29-00 

07-00 

08-00 

08-01 

10-00 

11-00 

12-00 

13-00 

14-00 

15-00 

20-00 

50-01 

51-00 

52-00 

09-OOt 

16-OOt 

19-OOt 

21-OOt 



;.176 



+0.298 
-0.401 



7.002 1°;^34 
11.030 tl-lfs 



6.105 



f0.192 
-0.204 



o IOC +0.037 
o.ioo _o.048 

1.0D4 _o.o51 

9 croQ +0.073 
^.ooo _o.091 



5 250 +°-^-'° 
O.ZOU _o.l32 

9.201 

8.282 



0.399 
0.500 
1.247 
1.567 

8.618 _i;752 

10.667 tlitl 



7.977 



+0.501 
-0.649 



o 091 +0.334 
o.ozi _o 473 

11.128 1°;1S 
11.501 i°:«^ 



S.221 



+0.152 
-0.322 



9.202 ±?;54| 

q n7Q +0.041 
O.UIJ _o.040 



f0.037 
-0.040 



4.126 

3 812 +0-127 

9 054 +°-*^^ 

o 010 +1.784 
o.o±o _o.969 

7 qcc +0.596 

10.971 t°;?|° 



o.»ou_Q gQ3 

n 041+0.439 

o.O'i±_Q 47Q 
2 224+°'^2^ 



3.742 



+0.197 
-0.204 



1 56q+°-2'*^ 

l.dU»_Q 236 



1.754 



+0.039 
-0.021 



1 Sfil +0-174 
l.»bl_(, 209 
o OOC.+0.157 

0.00<J_g 233 

9 584+°-*^°* 

».JO'i_Q 7gQ 

10 791+0.685 
lO./zig g42 

13.078l:i;734 

4 fjoc:+2.678 

'i.OOO_g ggg 

10.268l°;5?° 

Q oci +0.826 
O.ZO±_Q 269 

8.312+°;«J 
9 279+1-751 

».Z/y_(j ggg 

q i7n+1.269 
o.l/u_Q 728 

10.267t\ill 

0.678l°;l|° 

,+0.138 







-0.104 



9 77Q+0.299 
z./ lo_Q 247 

8.648l°2:°f 

10.529t°;f7°7 

5.82lt;i?^ 



Table 4. Frequency um, width Ajjj of the harmonic QPO Ljjj 
from the best fit to the PDS. 
1 2 sigma detection. 



Table 3. Frequency uu, width A// of the harmonic QPO Lji 
from the best fit to the PDS. 
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Obs ID 1^1/2 Ai/2 

10408-01-22-01 1.298 ^^'^2 l-lOSlg'.ng 

10408-01-22-02 1.255 t^'^s 0-96llo;ig4 

10408-01-23-00 2.026 Igogg 0.869t°;^°7 

10408-01-24-00 1.024 t'-Q^'^i^ 1.012J:|j-j^^ 

10408-01-25-00 0.532 tl'^i ^■^^'^-I'Xol 

10408-01-28-00 0.492 1^;[5|4 0.618lg;Jgg 

10408-01-29-00 0.882 t|j;[;^^ 0.919tg-j^^ 

10408-01-30-00 1.620 lo;i|g 2.229t°;2e6 



10408-01-31-00 1.556 _o 3(37 
20402-01-04-00 



0.241 ^ 904+0-440 

0.307 ^■^'-'^-0.430 
+0.065 1 qno+0-146 



20402-01-07-00 1.041 _o 03-. 
20402-01-08-00 1.406 



0.034 1 RflQ-l-0.098 



+0.046 1 fl«K+0-104 
-0.046 J^-°"'-'-0.098 
+0.055 1 Rfii +0.163 



0.030 n Rqn+0-150 



20402-01-08-01 1.277 _o ogj 

20402-01-09-00 1.346 tijljjgg 0.303l^;J|^ 

20402-01-10-00 1.440 Ig 029 

20402-01-11-00 1.363 Igol^ 1.872t°;^g^ 

20402-01-12-00 1.238 10027 '^■^'^^t^o^ist 

20402-01-13-00 1.575 l°;°2i l-^SOtoogg 

20402-01-14-00 1.567 tooj" l-956lo;i43 

20402-01-15-00 1.092 to'.ow 0-325loogg 

20402-01-18-00 1.584 ^JJ-fJ^ 0.684l|j-jf* 

20402-01-19-00 1.077 Igoj? 0-164J;|5:o63 

20402-01-20-00 1.592 llJUja 0-563l°;ogg 

20402-01-21-00 1.862 l°;5g 0-806J:ol57 

20402-01-49-00 1.127 tl'^°^\ "^■^^'^'^I'^^A 

20402-01-49-01 1.316 to;og4 LlOSlo'igy 

20402-01-50-01 0.510 tgloji 0-370lo;"| 

20402-01-52-00 0.744 t^Q^l 0.442l!^-22i 

Table 5. Frequency 1^1/2 and width ^1/2 °f the Lorcntzian ( 

ponent L1/2, from the best fit to the PDS. 
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ObsID 



i^h 



10408-01-22-01 


3.647 


10408-01-22-02 


3.398 


10408-01-23-00 


4.101 


10408-01-24-00 


2.647 


10408-01-25-00 


1.461 


10408-01-27-00 


0.571 


10408-01-28-00 


1.257 


10408-01-29-00 


2.595 


10408-01-30-00 


4.288 


10408-01-31-00 


4.831 


20402-01-04-00 


2.516 


20402-01-07-00 


4.606 


20402-01-08-00 


4.533 


20402-01-08-01 


4.776 


20402-01-09-00 


2.121 


20402-01-10-00 


1.884 


20402-01-11-00 


3.806 


20402-01-12-00 


2.842 


20402-01-13-00 


4.646 


20402-01-14-00 


4.082 


20402-01-15-00 


1.089 


20402-01-16-00 


3.612 


20402-01-18-00 


3.173 


20402-01-19-00 


1.432 


20402-01-20-00 


3.026 


20402-01-21-00 


3.533 


20402-01-49-00 


3.815 


20402-01-49-01 


4.033 


20402-01-50-01 


1.209 


20402-01-51-00 


1.662 


20402-01-52-00 


1.647 



+0.122 
-0.169 
+0.121 
-0.162 
+0.025 
-0.025 
+0.092 
-0.095 
+0.030 
-0.066 
+0.019 
-0.025 
+0.045 
-0.081 
+0.075 
-0.111 
+2.013 
-2.607 
+0.237 
-0.263 
+0.076 
-0.109 
+0.065 
-0.080 
+0.163 
-0.173 
+0.095 
-0.173 
+0.224 
-0.602 
+0.404 
-0.310 
+0.114 
-0.447 
+0.172 
-0.140 
+0.109 
-0.168 
+0.186 
-0.150 
+0.574 
-0.105 
+0.340 
-0.894 
+0.130 
-0.515 
+0.215 
-0.213 
+0.082 
-0.183 
+0.860 
-0.435 
+0.190 
-0.231 
+0.138 
-0.168 
+0.209 
-0.268 
+0.067 
-0.063 
+0.056 
-0.058 



1.569_Q 278 

1 71C-+0.333 
l.(10_Q 275 



0.932 

1.151 



,+0.035 
-0.037 
+0.120 
-0.176 



U.DZ(_g Qgg 

0.686t°;r7 
0.574t°:J«7l 

1 400+°-^™ 
1.1UU_Q 240 

19 979+0013 
ly.a I »_QQ52 

c 077+0. 242 
'-'■•^' ' -0.369 

0.959l°;188 

l.ZDl_o 253 

1 S9/1+0-241 
l-»^4_03g4 

1 052+0-549 
l.UO^_0 263 
o Q/lf;+l-230 

0.0'iU_Q gg2 

6 567+0-442 
0-00'_1.776 
, 779+0.917 
^■' '^-0.328 
+0.709 
-0.357 
,+0.662 
'-0.265 



3.020 
2 



1 SSQ+0-240 
i.8»y_0.188 

5 042+0-''35 

q 071+1.698 

•^■'" ^-1.003 

3 314+2-131 



4.392 



+0.617 
-0.659 



q 9Qc:+0.840 
o.z»o_Q 423 

4.1Oy_0.810 
4 170+0.182 
^-l '^-0.197 

4 061+°-^^° 

1.UDl_o j^gg 

2 679+0-11* 
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Table 6. Frequency j/^ and width A^ of the Lorentzian compo- 
nent Lft, from the best fit to the PDS. 
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